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ABSTRACT 

We explain the large scale correlations in radio polarization in terms of the cor¬ 
relations of galaxy cluster/supercluster magnetic field. Assuming that the po¬ 
larization correlations closely follow the spatial correlations of the background 
magnetic field we recover the magnetic field spectral index as —2.74 ± 0.04. 
This remarkably agrees with cluster magnetic field spectral index obtained in 
cosmological magneto-hydrodynamic simulations. We discuss possible physi¬ 
cal scenarios in which the observed polarization alignment is plausible. 
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1 INTRODUCTION 


The polarization directions from distant radio sources have been observed to show an align- 

ment over a distance scale of order 100 Mpc( 

Tiwari & Jain 

2013 

). This alignment is seen in 

the significantly polarized sources in JVAS/CLASS data ( 

Jackson et al. 

2007 

) with polar- 


ized flux greater than 1 mJy. Such a global alignment of polarization angles is unexpected 
but not in conflict with any fundamental principle. The distance scale 100 Mpc corresponds 
to the scale of galaxy superclusters and at such distances it is not unreasonable that the 
galaxies may show some correlation with one another. A similar dist ance scale also emerges 


in the study of galax y correlations using Sloan Digital Sky Survey flEisenstein et ah 


2005 


Anderson et al. 


20141) and agrees with predictions of the Big Bang cosmological model. How¬ 


ever a precise physical mechanism which might lead to an alignment of polarizations is so 
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far not available in the literature. Polarization measurements are affected by several instru¬ 
mental and observational biases which may also lead to the observed signal. We discuss 
some of these in section |3l An import ant observation is that the alignment is absent for 
low polarizations flTiwari fc JainI 1201311 . This provides some restriction on any explanation 
in terms of instrumental bias. Nevertheless, the issue of contribution due to bias requires 


a more detailed study. Furt hermore the sf 


gnihca nce of the alignment effect is found to be 


approximately three sigmas flTiwari fc Jainll2013[l . As we shall show in the present paper, the 


signihcance is further reduced if we take the jackknife errors into account. Hence the effect 
may also be simply due to a statistical fluctuation. On the other hand, assuming that the 
observations represent a real effect, it is interesting to think of possible physical processes 
producing such large scale global radio polarization alignments. 

In this paper we propose a physical process which may potentially explain the observed 
alignment. The model is based on two major assumptions: 

1. The galaxy jet axis is correlated with the cluster magnetic held. 

2. The integrated radio polarization is correlated with the jet axis. 


is iustihed bv several o 

3serva 

tions ( 

Gabi 

zda et al. 

1994; 

Lister & Smith 

2000; 

Pollack et al. 

2Q02; 

Helmboldt et al. 

2007; 

Joshi et al. 

2007 

) which indicate that the integrated polariza- 


tion from such sources is predominantly perpendicular or, less frequently, aligned with jet 
axis. Furthermore, the galaxies are known to be statistically aligned over large distance 
scal es, althoug h a pr o per physical understan ding of this phenomenon is so far lacking (see 


e.g. 


Kirk et ah 


(120151) : iKiessling et al.l (120151) . and references therein). The cluster magnetic 


held strength, as observed in cosmological magneto-hydrodynamic simulations , closely fol 


lows the cluster matter density prohle outside the core region of the cluster (iDolag et al. 


20021 ). The power sp ectrum of the mag netic held can be approximated by a power law with 
an exponent ~ —2.7 (lDolag et al.ll2002l) . Given that the galaxies show some alignment over 


large distances and the fact that cluster magnetic held is correlated with the matter density 
prohle may provide some motivation for our hrst assumption. The presumed magnetic held 
is expected to show some large scale correlations in real space. This is discussed in more 
detail in Section O Since, by our assumptions, the integrated radio polarizations are corre¬ 
lated with the background magnetic held, we expect the polarizations of diherent galaxies 
to be aligned with one another over the cluster or supercluster distance scale. 






























































Extracting Spectral Index of Intergalactic Magnetic Field from Radio Polarizations 3 
The JVAS/CLASS data sources are core-dominated flat spectrum radio source and are 
predominantly quasars and BL Lacs. The integrated polarization from these sources is a few 
percent. Here we consider only the significantly polarized sources for which the polarization 
is greater than 1 mJy. Based on the assumptions stated above, our physical mechanism 
implies that the integrated radio polarizations are correlated with the cluster magnetic 
held. Hence the observed alignment of the radio polarizations contain information about 
correlations of the magnetic held. We use this relationship in order to extract the spectral 
i ndex of the cluste r magnetic held, whose power spectrum is assumed t o follow a power law 


(IDolag et al 


2003 


200211. We point out that the primordial magnetic held (ISubramanian et al 


Seshadri fc Subramanian 


2005, 


2009 


Jedamzik et al 


1998; 


Subramanian Sz Barrow 


19981 ) is also expected to show a power law behaviour, however, in this case the spectral 


index is expec ted to be very di fferent in comparison to the expectation for the cluster 
magnetic held flDolag et al.ll2002[) . We simulate the cluster magnetic held for some assumed 
value of the spectral index. The correlations of the magnetic held directions at different 
spatial positions are assumed to be directly related to the corresponding correlations of 
the polarization angles, i.e. the latter provides an unbiased estimate of the magnetic held 
correlations. We study these correlations by defining a statistic Sn or 5^, as discussed 
in Section 01 By making a fit to the data statistic we extract the spectral index of the 
cluster magnetic held. The polarization data is likely to have large scatter and we use the 
jackknife estimate of errors. The jackknife errors are found to be large in compa rison to 


the squared variance of alignment statistics of shuffled PAs (ITiwari fc JainI 120131 ). In our 


earlier determination of the significance of alignment we had used the latter procedure 
flTiwari fc Jainll2013[) . If we instead use the jackknife errors we And that the significance of 
alignment is reduced. However, even with jackknife errors, the alignment signal is found to be 
good enough to sharply constrain the spectral index of magnetic held assuming the model 
presented above. We hope that the situation will improve with future Square Kilometre 
Array (SKA) observations. 

The paper is organized as follows. In Section [2] we describe the magnetic held model and 
explain its correlations in real space. We also explain our numerical procedure to generate 
a full 3D realization of magnetic held for a particular set of parameters. In Section [3] we 
give details of the JVAS/CLASS data and discuss the observed alignments. In Section 0] 
we review different statistical measures of alignment used in this paper. We describe our 
procedure in Section |5l We present our results in Section |6] and conclude in Section 0 
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2 CORRELATIONS IN INTER-GALACTIC MAGNETIC FIELD 


Magnetic field has been observed at all scales in inter-galactic medium. However the origin 
of ob s erved magnetic held is unknown and most likely to be primo r dial fISubramanian et ah 


2003 


Seshadri fc Subramanian 


2005, 


2009 


Jedamzik et ah 


1998; 


Subramanian fc Barrow 


19981 ). In our analysis we need to simulate the intergalactic magnetic held. On the cluster 


scale, cosmological magneto-hydrodynamic simulations suggest that the power spectrum of 
the m agnetic held can be modelled as a power law with an exponent ~ —2.7 flDolag et ah 


2 OO 2 I) . However on larger distance scales we expect that this exponent may be diherent. 


Hence we expect that a simple power law may not be valid at all distance scales. Here 
we assume a simple power law form of the power spectrum at all scales with an exponent 
corresponding to the cluster magnetic held. This will correctly reproduce the magnetic held 
correlations on the cluster scale, which is the only scale of interest in our analysis. It will fail 
at larger distances which are not of interest in the present work. Hence this failure cannot 
ahect our results. 

Let hi{k) represent the magnetic held in Fourier space. We can express its two point 
correlations as. 


{m)h(9))= ( 1 ) 
where k = \k\ and Rj is the projection operator given as, 

Pi = (% - . ( 2 ) 

The real space magnetic held can be written as, 

= (3) 

k 

where V is the volume. We assume a power law dependence of the spectral function M{k) 


M{k) = Ak'^^, 


(4) 


with the spectral index ub > —3. The magnetic held is assumed to be statistically uncorre¬ 
lated in fc-space. However the held is correlated in real space and the nature of correlation 
is controlled by the index ns- We can write the real space correlation of the held as the 
Fourier transform of equation ([1]). We obtain, 

(R(r+ f')R,(F)) = ^j dV-^'P,,{k)M{k)W\krG), (5) 



































Extracting Spectral Index of Intergalactic Magnetic Field from Radio Polarizations 5 
where is the ’’galactic” scale taken as 1 Mpc and W is a. window function of the form, 


W{x) = 


1 X < 1 
0 X > 1 


( 6 ) 


This window function hxes the value of the scale re- The magnetic held B{r) is assumed 
to be uniform over the scale r^ = = 1 Mpc. In equation ([5]) we have also taken the 

continuum limit and replaced V I constant A in equation (|4]) is equal to 

■ It is hxed by demanding {Bi{r)Bi{r)) = Bq, where Bq is the intergalactic 

kfj 

magnetic held averaged over the dis tance scale of rr; = 1 Mpc. We expect that for the 


case of the primordial held, Bq ~ nG (jSeshadri fc Subramanian 


De Angelis et ah 


2009; 


Yamazaki et ah 


2010 : 


20081) ■ although this value plays no role in our analysis. Furthermore we 


add that the radio sources considered in this work are separated from one another by a 
mean distance of tens of Mpc (Fig. [2]) and, hence, averaging magnetic held over 1 Mpc is 
reasonable. 

It is also appropriate to demand a large scale cut-oh for the correlations in equation (j5]) . 
Here we assume that this cutoh is sufficiently large (r^ax >3Gpc) so that we can simply set 
^max to be oc. Heuce, we set the lower limit of integration in equation ([5]) as kmin = ''"mL = 0- 
We emphasize that such large scale correlations are expected within the Big-Bang cosmology. 
The perturbations at the time of inhation have wavelengths larger than the comoving size 
of the current observable universe and hence the primordial magnetic held correlations may 
also exist over the horizon scale. 

We num erically generate inte rgalactic magnetic held in hD space using the procedure 


described in 


Agarwal et al.l (120121 ) for diherent values of the spectral index hb- We consider 


discretized space consisting of a large number of cells or domains of equal size. The magnetic 
held is assumed to be uniform in each domain. We hrst generate the magnetic held in k- 
space, which is straightforward since the corresponding held is uncorrelated. It is convenient 
to use polar coordinates (/c, 6, 0) in fc-space. The projection operator Pij{k) ensures that 
the component of magnetic held along k, bk is zero. The remaining two orthogonal compo¬ 
nent bg and b^ are uncorrelated and therefore we generate these by assuming the Gaussian 
distribution. 


f{bg{k),b^{k)) = N exp 


_ I bj{k) + bl{k) \ 

V 2M(fc) J 


(7) 


where N is the normalization. The distribution in equation ([7]) represents an uncorrelated 
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magnetic field in fc-space. Next, we do a Fourier transform to obtain magnetic field in real 
space. 


3 DATA 

We use the catalogue produced by 


Jackson et al. 


(120071) . The catalogue contains 12743 core 


dominated flat spectrum radio sources and provides their angular coordinates and the Stokes 
/, Q and U parameters. The observable required for the alignment study is the polarization 
a ngle (PA). A de t ailed description of the catalogue including the calibration methods is given 
(120071 ). We hnd that there are 4400 sources in catalogue with polarized flux 


m 


Jackson et al. 


greater than 1 mJy. JVAS/CLASS sources are presumably located at very large distances 
{z ~ 1), although the exact redshift of each source is unknown. A signal of alignment 
of radio polarization angles (PAs) has been found in signihcan tly polarized sources , with 
polarized flux greater than 1 mJy, contained in this catalogue (iTiwari fc JainI 120131) . The 


(Tiwari & Jain 

2013) 

result in 

Joshi et al. 


also indicates alignment in data (IShurtlefB 120141) . Furt hermore alignments have 


Deen found 


even at larger distances for QSOs in this data sample (jPelgrims fc Hutsemekers 


20071: 


le data may 


3e aff ected by some instrumental and observational biases 


Jackson et al 


2015h. 


Joshi et al. 


20071) . One possibility is the error in the removal of residu al instrumental 


polari zation. This may artihcially generate large scale alignments. However 


Tiwari fc Jain 


(1201 3l ) argue that this must dominate for sources with low polarizations which do not show 
any signal of alignment. Hence it is not possible to attribute the observed alignment to this 
bias. Another possibility is that sources in a small neighbourhood are observed together 
within a particular observational run. It is possible that this could generate alignment in 
sources within small angular separations, as observed in iTiwari fc JainI (120131) . We cannot 
rule out such a possibility. This issue is best addressed by future more rehned observations. 

In our analysis we need redshift distribution of JVAS/CLASS sources in order to model 
the statistical alignments of radio polarization. Since the redshifts of many sources are un¬ 
known, we adopt the following hybrid redshift model. We employ NASA/IPAC EXTRA- 
GALACTIC DATABASE (NED) —‘Retrieve Data for Near- Object/Position List 3 tool 


^ https://ned.ipac.caltech.edu/forms/rmd.html 



















































Extracting Spectral Index of Intergalaetic Magnetic Field from Radio Polarizations 7 
and cross match all (above 1 mJy) 4400 JVAS/CLASS sources with NED objects and hnd 
nearby objects within radius 0.1 and 0.5 arcsec. We identify 1783 sources with a search 
criteria of 0.1 arcsec and an additional 138 if the larger radius of 0.5 arcsec is used. This 
leads to a total of 1921 sources with redshifts in NED catalogue. Most of these (1389) are 
quasars, a few (201) are listed as galaxies and only 11 are identihed as radio sources. The 
remaining 320 sources are other types or unknown type. We show the redshift distribution 
of all these different classes of sources in Fig. [1] Even so we do not hnd the redshift for 2479 
(out of 4400) sources and we adopt the radio source redshift prohle for these. This is justihed 
because the redshift for radio sources is largely unknown and thus most likely the unknown 
redshift sources are radio only. Anyway, in our analysis we do not need precise redshifts, 
we only need redshift distribution to glean out the magnetic held statistical correlations. 
The radio redshift distribution is largely unknown and we only have a few small area deep 
survey observations to determine the redshift number density. We rely on Comb ined EIS- 


NVSS Survey of Radio Sources 


Hercules (jWaddington et al. 


( CENS ORS! flRest et al. 


2000 


2003 : 


Rigby et al.lboiD and the 


200lh observed redshift number distribution and model 


the radial nu mber density of remaining 2479 sources. This is 


as followed in 


Nusser fc Tiwaril f 2015 h iTiwari fc Nusserl ( 2015 ) and presumably the best we 


he same redshift distribution 


can assume for JVAS/ CLASS radio sources. W e show the CENSORS and Hercules redshift 
number distribution ht (iNnsser fc Tiwarill2015l ) as ‘ht’ in Fig. [T], a sample redshift histogram, 
input to simulation, is also shown in the same hgure. 

We remove 187 sources from these 4400 sources as their angular positions coincide with 
other sources. We point out that the sources lie dominantly in the Northern hemisphere. 
Furthermore there are very few sources along the galactic plane. We use only this sample of 
4213 sources for our analysis. 

The alignment statistic we use is dehned for the number of nearest neighbours, n^, or 
equivalently the angular separation A6 between sources. For any chosen source k we order 
the remaining sources in terms of increasing angular separation from the source k. The 
closest Uy sources, excluding the source k itself, forms the required nearest neighbour set. 
The minimum value of is clearly 1 and the maximum value we explore is equal to 15 for 
reasons given later in section [6l We can associate a mean A6 with Uy nearest neighbours by 
determining the A6 corresponding to the nearest neighbour set of each source and taking the 
average over all sources. The mean value of A6 increases monotonically with Uy. Hence Uy 
also provides a measure of the angular separation. In order to relate Uy or A9 to separation 
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Figure 1. The redshift distributions of different classes of JVAS/CLASS objects. A total of 1921 sources out of 4400 are 
retrieved from the NED database. The filled (or dashed line) histogram corresponds to NED object match within 0.1 arcsec, 
whereas the empty (or solid line) corresponds to NED object match within 0.5 arcsec. Assuming that the remaining 2479 
sources follow radio source redshift number density profile, we have shown one sample input to our simulations. 


Distance (Mpc) 

60 80 100 



Figure 2. The average angular separation for a fixed number of nearest neighbours {riv) and vice versa. Corresponding physical 
distance scale (Mpc) of the source assuming that it is located at 2 = 1 is also given. 


distance we assume that sources are located at mean redshift approximately equal to 1. We 
need to make this assumption since we do not have redshifts for most of sources, indeed 
we only have the probability distribution. We have the redshifts for 1921 sources, most of 
them are quasars and are foun d to have larger distances scale correlations in 2D and 3/1 
fjPelgrims fc Hutsemekersll2015h . We show the average angular separation and the physical 
distance scale as a function of in Fig. [2l The physical distance scale L of a distant cluster 
is related to its angular size A6 by the formula L = JaAO, where Ja is the angular diameter 
istance, compu ted at redshift 2 ; = 1 using the s t anda rd Lambda Cold Dark Matter model 


(IWeinbergi 120081 ) . We clarify that iTiwari fc JainI (120131 ) used the comoving distance instead 


of the angular diameter distance. This leads to a change in our estimate of the distance scale 
of alignment from 150 Mpc to 100 Mpc. 
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4 CORRELATION STATISTIC 


In order to quantify the correlations between polarizations at different positions we define 
two statistics, Sn and S'^. Consider the nearest neighbours of a source located at site 
k. Let ipi be the PA of the source at the site within the nearest neighbour set. The 
dispersion of PAs in the neighbourhood of this site may be characterized by the measure, 


dk 




cos[2('0i + Ai^k) - 2t^fc)]. 


( 8 ) 




Here the factor Ai^k arises since the polarizations at two different points, labelled as i and 
k, on the celestia l sphere have to be correlated after making a parallel transport from i ^ k 
(jjain et al.ll2004[) along the geodesic which connects the two positions. The parameter in 
this equation is the PA at the site k. The measure dk provides an estimate of the dispersion. 
We point o ut that a large value of dk implies lo w dispersion and vice versa. The statistic is 


dehned as (iHutsemekers 


1998 


Jain et ah 


200J, 


SD = -^dk, 


(9) 


k=l 


where Ug is the total number of data samples. A strong alignment between polarization 
vectors implie s a large value of Sn - We point out that the notation used in this work is 


different from iTiwari fc JainI (l2013f) . We have also dropped self correlations {i = k) while 
calculating dk in equation (|8]). This changes the magnitude of Sd signihcantly for small 
but does not change the alignment signihcance results. However, as already mentioned above, 
in this paper we are considering the jackknife errors which are signihcantly larger than the 
varianc e oi Sn for the case of the random polarization samples, considered in our previous 
paper flTiwari fc Jairul2ni3h . This does lead to an appreciable change in the signihcance of 
alignment. 

Alternatively we dehne statistic S'j^ where we measure the dispersion over a hxed angular 
separation {A9) rather than measuring it over a hx number of nearest neighbours (u„). We 
point out that these two statistic Sr, and S'j^ will give same result for a spatially uniform 
source distribution since for a circle of a given radius the nearest number of sources will be 
same everywhere. Hence, Sd can be translated into S'r) naively. However, the data sample 
is expected to show some deviation from spatial uniformity and hence we expect small 
diherences in the two measures, Sd and S'r). 
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5 PROCEDURE 


We first determine the statistic Sd ( see equation ([9])) for the observed linear PAs. The same 
procedure is repeated for the theoretical PAs, which are assumed to be aligned perpendicular 
to the background magnetic held. The magnetic held is generated by simulations as explained 
in Section [2] and we glean out the galaxies according to redshifts details as given in Fig. [T] 
For sources where we only have probability distribution of number density, we sample over 
100 random outputs of redshifts generated from the probability distribution shown in Fig. 
m Finally, we calculate the 2D correlation statistics Sd and S'j^. The resulting theoretical 
values oi Sd are computed for a range of values of the spectral index ub- The best ht value 
of is obtained by making a fit fo fhe observed data. 

In order to compute ^6 need an estimate of the error in the statistic Sd- We resort 
to jackknife errors as we do not know the exact errors in polarization measurements. For 
computing the jackknife errors, we resample the data by eliminating the source and 
calculate the correlation statistics SD^i)- We repeat this process for all sources and calculate 
the SD{i) for z = 1 to 4213. We call the full sample statistics as Sd and the jackknife error 
in its estimation is given as 

\ N 

(SSnf = (10) 

where N is the total number of data point, which is 4213 in our case. We add that the 
jackknife errors do not include systematic errors. The distribution of jackknife sampled Sd 
and a Gaussian £t for = 10 is shown in Fig. [31 Note that the mean of jackknife sampled 
statistics is almost the same as the full sample Sd and so the standard deviation in Fig. [3] 
is roughly equal to 6Sd/ \/(A^ — !)• The SSd for = 10 as estimated from equation (IT0|) 
is found to be 0.007716. We similarly calculate the jackknife errors for all and angular 
scales A6 for statistics Sd and S'j^. In Fig. [Hand[5] the jackknife errors are shown in lowest 
panel. 


6 RESULTS 


We simulate the 3D magnetic field for ub values in range —2.20 to —2.98, considering the 
power law spectrum as discussed in Section [2l As described in Section [31 we have redshifts of 
19 21 sources and fo r the r emaining sources the redshift is generated assuming a fit obtained 
in [Nusser fc Tiwaril (1201511 . We glean out the source location magnetic held directions from 
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Figure 3. The distribution of jackknife sampled statistics Sd for Hv = 10. A fit to Gaussian is also shown. Note that the full 
sample statistics Sd = 0.01106 is almost equal to the jackknife sampled statistics Sd distribution mean (0.01106) and so the 
standard deviation in figure is roughly the SSd/\/{N — 1). 



Hv 


Figure 4. The statistic Sd for data with jackknife error bars along with Sd corresponding to randomly generated PAs (open 
circles) are shown in the top panel. The alignment significance considering jackknife errors is plotted in middle panel. For 
clarity we have plotted jackknife errors SSd for each riy in bottom panel. 


simulated magnetic field and calculate the model statistic Sd- We average over 100 real¬ 
ization to sample over fit generated random redshift positions. The data statistic Sd along 
with alignment signihcance and jackknife errors is shown in Fig. 01 The model Sd values for 
different ub are htted with observed data Sd- The resulting is dehned as 

’Fb (observed) — Sd (simulated) \ ^ 

^ ) ’ 


= 


15 


n„ = l 


( 11 ) 


where 6Sd is the jackknife estimate of error in S'D(observed) (see Fig. 0] and [5]). We similarly 
calculate ^d- The maximum value of is set equal to 15. This choice is made 




















12 Prahhakar Tiwari and Pankaj Jain 



Figure 5. The statistics S'^. Data with jackknife error bars along with corresponding to random PAs (open circles) are 
shown in the top panel. The alignment significance considering jackknife errors is plotted in middle panel. The jackknife errors 
are in bottom panel. 


since the significance of alignment beyond = 15 is within 1-a. With one parameter the 
number of degrees of freedom (DOF) is equal to 14. We present the /DOF values verses 
spectral index ns for statistic Sn in Fig. O The minimum value of /DOF is found to be 
approximately 0.25 at ub = —2.76. The data and best fitted simulated Sb comparison for 
spectral index ub = —2.76 is given in Fig. [71 Including one sigma error, the extracted value 
of Ub is found to be, ub = —2.76 ± 0.04. We point out that the relatively low value of 
indicates that the jackknife errors are large. Nevertheless, even with such large errors we are 
able to clearly resolve the magnetic field spectral index as seen in Fig. [6l 

The results for the case of the alternate statistic S'^ are also shown in Fig. O In this case 
we set the maximum value of the fixed angular distance A6' such as to include 15 nearest 
neighbour as an average. This corresponds to A0 = 4.4°. The full alignment results for 
statistic S'b are shown in Fig. [5l The angular distance corresponding to a mean value of 
is somewhat higher that the average distance for a fixed value of n^. We have shown this 
difference in Fig. [2l Despite this difference, we obtain almost similar results for statistics 
S'b- The best fit value of the spectral index is found to be ub = —2.74 ± 0.04 (Fig. [6]). The 
slight difference in the extracted values oi ub for the two statistics reflects the spatial non- 
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uniformity of data. The value extracted using the statistic S'jj may be more reliable since it 
includes sources within a hxed distance from a particular source rather than including a hxed 
number of nearest neighbours. The /DOF is slightly lower for S'j^ but the data points at 
low angular separation (A0 < 2°) show large fluctuation due to spatial non-uniformity of 
data. Nevertheless, the results from So and S'j^ agree within errors. It is interesting that the 


extracted value of no is in good agreement with t hat obtainec 


simulations which suggest a spectral index of —2.7 flDolag et al 


by m agneto-hydrodynamic 


2 OO 2 I) . The close agreement 


may be fortuitous but may be tested by future more refined data. In any case, theoretically, 
we do not expect a perfect agreement between these two indices. The polarization index {ns), 
extracted from observations, only acts as a tracer of cluster magnetic field index If we 
assume that the correlations in polarization are induced by those in the cluster magnetic 
field, we expect that the maximum level of alignment in the polarization orientations would 
be equal to those of the background magnetic held. The level of alignment, i.e. the value of 
the statistic So or S'o, increases with \nB\- Hence we expect that \nB\ ^ Wb\- 

We also point out that since the best ht value of the spectral index is close to the 
theoretical expectations, we could not have expected a higher signihcance of alignment in 
this data set. If we compute the signihcance using the variance of alignment statistics of 
shuffled PAs flTiwari fc JainI 120131 ). it is found to be about 3 sigmas. However if we include 
the jackknife errors, the signihcance is at best about 1.5 sigmas, as shown in Figs. |4]and|5l 
In order to increase it to 3 sigmas we need to reduce the errors by a factor of two which 
requires four times more data. Hence a clear test of our proposal can be made by acquiring 
a suitably enhanced data set. 


7 CONCLUSION AND DISCUSSION 

We have presented a possible model to explain the large scale radio polarization correlations. 
The model is based on two main assumptions that the galaxy jets are aligned with cluster 
magnetic held and that the jet orientations approximately mark the radio polarization an¬ 
gles. We argue that the second assumption is well supported by data and the cosmological 
magneto-hydrodynamic simulations provide some support for the hrst assumption. We also 
hnd that our extracted value of the spectral index {hb = —2.74 ± 0.04) is in good agree¬ 
ment with the value —2.7 obtained by magneto-hydrodynamic simulations. This is rather 
encouraging and provides additional support to our proposal. However given the inherent 
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Figure 6. The jDOF as a function of the spectral index ng for the statistics and 5^. 

uncertainties in the polarization data, we cannot claim this to be dehnitive evidence for our 
model, which requires further testing with more rehned data. We conclude that the observed 
alignment in the JVAS/CLASS data can be successfully explained in terms of magnetic held 
correlations. The model can be further applied to other data sets and tested. 
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